High Resolution X-ray Spectroscopy of the Post-T Tauri Star 

PZ Tel 



C. ArgirofR^, J.J. Drake^, A. Maggio'^, G. Peres^, S. Sciortino^, and F.R. Harnden^ 



O 

o 

(N 



> 

O 
cn 
o 

o 

:^ 

o 
U 



X 



ABSTRACT 

We present an analysis of the Chandra High Energy Transmission Grating Spectrometer 
observation of the rapidly rotating (Prot = 0.94 d) post T Tauri (~ 20 Myr old) star PZ Telescopii, 
in the Tucana association. Using two different methods we have derived the coronal emission 
measure distribution, em{T), and chemical abundances. The em{T) peaks at logT = 6.9 and 
exhibits a significant emission measure at temperatures logT > 7. The coronal abundances are 
generally ~ 0.5 times the solar photospheric values that are presumed fairly representative of 
the composition of the underlying star. A minimum in abundance is seen at a first ionization 
potential (FIP) of 7-8 eV, with evidence for higher abundances at both lower and higher FIP, 
similar to patterns seen in other active stars. From an analysis of the He-like triplet of Mg XI 
we have estimated electron densities of ^ 10^^ — 10^^ cm^"^. All the coronal properties found 
for PZ Tel are much more similar to those of AB Dor, which is slightly older than PZ Tel, than 
to those of the younger T Tauri star TW Hya. These results support earlier conclusions that 
the soft X-ray emission of TW Hya is likely dominated by accretion activity rather than by a 
magnetically-heated corona. Our results also suggest that the coronae of pre-main sequence stars 
rapidly become similar to those of older active main-sequence stars soon after the accretion stage 
has ended. 

Subject headings: stars: abundances — stars: coronae — stars: individual (PZ Telescopii) — stars: 
pre-main-sequence — techniques: spectroscopic — X-rays: stars 



1. INTRODUCTION 

One of the primary characteristics of low mass 
pre-main sequence (PMS) stars is their intense X- 
ray activity. This X-ray emission therefore rep- 
resents an important means for investigating the 
properties and evolution of young stellar objects. 
X-ray activity is present during the evolution of 
PMS stars both in the initial evolutionary stages 
of a Classic T Tauri Star (CTTS, Class I and 
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II sources), during which the star has an accre- 
tion disk that surrounds it, and in the subse- 
quent Weak-Line T Tauri Star (WTTS, Class HI 
sources) phase in which the star has no accret- 
ing material and is approaching the zero-age main 
sequence (ZAMS, Feigelson & Montmerle 1999). 
How the presence of accreting material influences 
the X-ray emission of PMS stars, and for how long, 
remain questions of debate. 

The high resolution X-ray spectra now available 
with the Chandra and XMM-Newton satellites of- 
fer the possibility to perform detailed studies of 
stellar coronae because key emission lines diagnos- 
tics can now be resolved. These diagnostics can be 
used to derive elemental abundances, temperature 
and density structure of the emitting plasmas. 

It is worth noting that very few young stars in 
star forming regions or associations are sufficiently 
X-ray bright to allow high resolution X-ray spec- 
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troscopy with current instrumentation. In partic- 
ular, among the CTTSs, TW Hya is the best stud- 
ied case to date because it is the nearest (~ 56 pc) 
known CTTS. TW Hya shows spectral character- 
istics very different from those of young but oth- 
erwise normal active stars (Kastner et al. 2002; 
Stclzcr & Schmitt 2004): very low plasma tem- 
perature (logT ~ 6.5), high density (log A^e ~ 13), 
very low Fe abundance (Ape/^Fe© ~ 0.2). Kast- 
ner et al. (2002) and Stclzcr & Schmitt (2004) at- 
tributed these characteristics to an accretion shock 
rather than coronal activity. Without the benefit 
of high resolution spectra of similar stars for com- 
parison, however, the nature of the peculiarity of 
TW Hya remains uncertain. 

The situation with slightly more evolved stars 
is more clear. Studies of AB Dor, a young ac- 
tive star that has nearly arrived at the ZAMS and 
which has been observed in the past with several 
space-borne X-ray observatories, have shown that 
it is characterized by a hot corona (T ~ 10^ K) 
with plasma densities ranging from 6 x 10^" cm~^ 
at 2 X lO^'K, to 3 X lO^^pj^j-a lO^K, and a 
moderately low Fe abundance {Ape/ApgQ ~ 0.25, 
Mewe et al. 1996; Giidel et al. 2001; Sanz-Forcada, 
Maggio, & Micela 2003b; P. Testa, in preparation; 
D. Garcia- Alvarez, in preparation). In many re- 
spects, AB Dor can be considered the prototype 
of very active single stars. 

High resolution X-ray spectroscopy of other 
PMS stars with ages between few 10^ yr, typical 
of CTTSs like TW Hya, and 10* yr, the ZAMS 
of solar-type stars like AB Dor, is crucial for un- 
derstanding how the characteristics of stellar X- 
ray activity change during these early evolutionary 
phases. 

In this work, we present a Chandra High Energy 
Transmission Grating Spectrometer (HETGS) ob- 
servation of PZ Telescopii (HD 174429, HIP 92680). 
Classified as a KOV star by Houk (1978), PZ Tel 
was determined by Zuckerman & Webb (2000) to 
be a member of the Tucana association, a nearby 
star forming region about 45 pc away. From ob- 
servations with the ROSAT PSPC, Stclzcr & 
Neuhauser (2000) deduced a PZ Tel X-ray lu- 
minosity of Lx ~ (2.88 ±0.08) x 10^° erg 8"^ At a 
Hipparcos parallax distance of 49.7 ± 2.9 pc (Per- 
ryman et al. 1997), PZ Tel is a single star with 
a rotational period of 0.94 d (Coates et al. 1980; 
Innis, Coates, & Thompson 1984; Innis, Thomp- 



son, & Coates 1986). Favata et al. (1998) have 
estimated a mass of 1.1 Mq for PZ Tel and an age 
of approximately 20Myr. Its youth is also con- 
firmed by prominc;nt Ha emission and a relatively 
undepleted Li abundance (Soderblom, King, & 
Henry 1998). Further evidence of its PMS status 
has been pointed out by Barnes ct al. (2000) who 
deduced from vsiui and Prot that the minimum 
radius of PZ Tel is R sini ~ 1.277?©: this value 
is larger than the radius of a main sequence star 
with the same mass as PZ Tel. 

The analysis of Chandra HETGS high reso- 
lution spectra of PZ Tel offers us the opportu- 
nity to study the coronal properties of a single 
star which has dissipated its accretion disk and 
is approaching the ZAMS. This study also allows 
us to compare the coronal properties of PZ Tel 
with those of both younger and older stars, pro- 
viding a glimpse of the evolution of stellar coro- 
nae through the PMS phase. In particular, the 
CTTSs have disks from which active accretion is 
still taking place. These stars could have quite dif- 
ferent magnetospheric geometries, possibly involv- 
ing magnetic connections between star and disk 
(e.g. Montmerle 2002). In the case of coeval and 
older stars, it appears that coronal activity of PMS 
stars without an accreting disk can be explained 
on the same basis as that of main-sequence stars 
(Flaccomio, Micela, & Sciortino 2003) — stellar 
rotation and convection. In this domain, PZ Tel 
provides a new window on phenomena such as the 
chemical fractionation of elements that is seen to 
occur in coronae over a wide range of activity level. 
In the solar case, elements with low first ionization 
potential (FIP) such as Mg, Fe and Si are seen to 
be enhanced relative to elements with high FIP, 
such as O, Ne and Ar (e.g. Fcldman 1992). In 
more active stars, the situation appears somewhat 
reversed, with elements such as Ne appearing en- 
hanced relative to those with lower FIP (e.g. Drake 
2002, and references therein). The case of the ac- 
tive but very young post-T Tauri stars remains 
unexplored at high spectral resolution. 

In § 2 and 3 we describe the Chandra observa- 
tion and the techniques used in its analysis. Sec- 
tion 4 presents the resulting coronal temperature 
structure and abundances. These are discussed in 
§ 5, in which we also present a comparison of the 
coronal properties of PZ Tel with those of AB Dor 
and TW Hya. 
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Fig. 1. — Light curve of PZ Tel obtained from the HEG and MEG spectra (excluding zero-order events) 
with bin size of 1000 s. 



2. OBSERVATION 

PZ Tel was observed with the Chandra HETGS 
on 2003 June 7, for 73.9 ks. The HETGS pro- 
vides both High Energy (HEG) and Medium En- 
ergy Grating (MEG) spectra. The spectral ranges 
available with these instruments are 1-18 A (HEG) 
and 1-25 A (MEG), with spectral resolutions A A 
of 0.01 A and 0.02 A respectively. We have used 
CIAO V3.0 to extract the positive and negative 
first order spectra for each grating and to compute 
the effective area of the HETG+ACIS-S combina- 
tions. In this latter step, we have also taken into 
account the attenuation resulting from the build- 
up of a contamination layer on the ACIS-S filter 
(e.g. Plucinsky ct al. 2003) using CIAO Version 1 
of the effective area contamination correction^. In 
Fig. 1 we show the light curve of photon arrival 
times of the HEG and MEG spectra, excluding 
zeroth order events that are subject to significant 
pile up. Figure 1 shows the X-ray flux to be char- 
acterized by variability on time scales of a few to 
tens of ks, with two flare-like events occurring at 
about 20 ks and 55 ks after the start of the obser- 
vation. 

In Fig. 2 we illustrate the first order MEG 
spectrum; labels identify the strongest emission 
lines. The X-ray luminosity obtained from the 
first order dispersed events in the interval 6-20 A is 

http: / / cxc.harvard.edu/cicio/threads / aciscontam / 



Lx ~ 2.2 X 10'^*'ergs ^, from which we derive 
log(Lx/iboi) = -3.2, having adopted Lboi = 
3.6 X 10'^^ ergs~^ (obtained from the my and B — V 
Hipparcos data and the bolometric correction of 
Flower 1996). 
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Fig. 2. — PZ Tel MEG smoothed spectrum (summed positive and negative orders) with bin size of 0.005 A. 



Table 1 
MEG Spectrum of PZ Tel. 
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Table 1 — Continued 
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0.19 


27.8 


1 


40b 




13.52 


FeXIX 


2p3 C^D) 3d 2s= 2p* 
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0.98 


± 
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Fo XVIII 


2p* (3p) 3d ^P^,2 2s^ 2p= 2p3/2 


6.90 


0.19 


± 


0.12 


23.4 




48b 




14.61 


Fo XVIII 


2p-» (3p) 3d ^P3/2 ^ 2s^ 2p= ^Pi/2 


6.90 












49a 


14.67 


14.67 


FoXIX 


2p3 (^D) 3s 2s^ 2p* ^P2 


6.90 


0.27 


± 


0.13 


23.7 


1 


49b 




14.67 


Fo XVIII 


2p^ CP) 3d ^D3,2 2s2 2p5 2Pi/2 


6.90 












50 


15.01 


15.02 


Fe XVII 


2p^ 3d ^Pi 2p^ ^So 


6.75 


1.85 


± 


0.23 


19.4 


1.2 


51 


15.08 


15.08 


FeXIX 


2s 2p* Cp) 3s ^P2 2s 2p^ ='P2 


6.90 


0.31 


± 


0.13 


20.0 


1 


52a 


15.18 


15.18 


O VIII 


4p ^P3/2 ^ Is ^5i/2 


6.50 


0.29 


± 


0.13 


19.5 


1 


52b 




15.18 


o VIII 


4p ^Pl/2 ^ Is '5i/2 


6.50 












52c 




15.20 


FeXIX 


2p3 CS) 3s ■'■Sa ^ 2s2 2p'' ■'P2 


6.90 












53 


15.26 


15.26 


Fe XVII 


2p^ 3d ^Di 2p^ ^So 


6.75 


0.44 




0.14 


18.1 


1 


54a 


16.00 


16.00 


Fe XVIII 


2p* (3P) 3s =P3/2 ^ 2s^ 2p= 2p3/2 


6.90 


1.00 


± 


0.17 


16.4 


1 


54b 




16.01 


O VIII 


3p ^P3/2 ^ 1« ^Si/2 


6.50 












54c 




16.01 


O VIII 


3p ^Pl/2 ^ Is ^Si/2 


6.50 












55 


16.08 


16.07 


Fe XVIII 


2p* (^'P) 3s *P5/2 ^ 2s^ 2p= ^P3/2 

2p=' CD) 3p ^P2 ^ 2s 2p^ =*P2 


6.90 


0.52 


± 


0.14 


16.1 


1 


56 


16.11 


16.12 


FeXIX 


6.90 


0.04 


± 


0.09 


16.0 




57 


16.16 


16.17 


Fe XVIII 


2s 2p5 (3p) 3s ^P3/2 2s 2p«' ^Si/2 


6.90 


0.25 


± 


0.12 


15.7 


1 


58 


16.77 


16.78 


Fe XVII 


2p^ 3s 3 Pi ^ 2p«' ^So 


6.70 


0.83 


± 


0.16 


12.4 


1 


59 


17.05 


17.05 


Fe XVII 


2p^ 3s ^Pi ^ 2p«' ^Sa 


6.70 


0.96 


± 


0.17 


11.7 


1 


60 


17.09 


17.10 


Fe XVII 


2p^ 3s ^P2 ^ 2p«' ^Sa 


6.70 


0.91 


± 


0.17 


11.1 


1 


61a 


18.97 


18.97 


O VIII 


2p ^P3/2 ^ 1« ^S'1/2 


6.50 


2.44 


± 


0.24 


7.2 


1,2 


61b 




18.97 


O VIII 


2p ^Pl/2 ^ 1« ^S'1/2 


6.50 












62 


21.60 


21.60 


O VII 


Is 2p ^Pi Is^ ^So 


6.30 


0.13 


± 


0.10 


3.3 


1,2 


63a 


24.77 


24.78 


N VII 


2p ^P3/2 ^ Is ^Si/2 


6.30 


0.12 


± 


0.10 


3.0 


1 


63b 




24.78 


N VII 


2p ^Pl/2 ^ Is 'S'i/2 


6.30 













^ Observed and predicted (CHIANTI database) wavelengths in A. In case of unresolved blends, identified by the same label number, we list the 
main components in order of increasing predicted wavelength. 

^ Temperature in K of maximum emissivity. 

^ Count rates in ctsks"^ with uncertainties at the 68% confidence level. In the cases of unresolved blends, identified by the same label number, 
we report only the total line count rates. 

MEG first order effective area in cm^, needed to convert the measured count rate to flux. 

*^ Lines selected to derive the emission measure distribution with approach 1 or 2. 
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3. ANALYSIS 

The MEG first order spectrum has more than 
1.6 X 10'* total counts and contains fines with good 
signal-to-noise ratios. The HEG spectrum con- 
sists of 5.5 X lO^cts, with few lines having suf- 
ficient counts to provide useful fluxes. We have 
therefore based our analysis on the MEG spectrum 
only. The MEG background is extremely low and 
contains an average of only 0.02 counts per ACIS 
pixel, and it was therefore ignored. 

The spectral analysis of PZ Tel has been per- 
formed using the PINTofALE VI. 5 software pack- 
age (Kashyap & Drake 2000). We have adopted 
the CHIANTI V4 emission line database (Young 
et al. 2003) and the Mazzotta et al. (1998) ioniza- 
tion equilibrium. 

Spectral line fluxes were obtained by fitting the 
counts histograms with a modified Lorentzian pro- 
file of the type 



/(A) - 



Ac 
AA 



with /3 = 2.4 plus piece- wise constant contin- 
uum level. For the continuum placement wc have 
adopted the following procedure. As a first step 
we have guessed an emission measure distribution 
em(T) and set of abundances in order to pre- 
dict a continuum level for measuring line fluxes. 
Prom these line fluxes we evaluated the em(T) 
and abundances for PZ Tel, and updated the pre- 
dicted continuum level. Considering this new pre- 
dicted continuum we have once again measured 
the line fluxes. Wc performed several iterations 
of the above procedure until convergence, i.e. un- 
til the continuum used for the line fitting agrees 
with the continuum predicted on the basis of the 
reconstructed em(T) and computed abundances. 

In Table 1 we report the strongest lines identi- 
fied in the PZ Tel spectrum, together with their 
measured MEG count rates. As stated in § 2, 
significant variability of the X-ray emission oc- 
curred during the observation, but the MEG signal 
to noise ratio did not allow us to perform time- 
resolved spectral analysis of the data. On the 
other hand, the integrated PZ Tel spectrum pro- 
vides a good representation of the average state of 
an active corona where a low-level flare-like activ- 
ity is present at all times. 



The hydrogen column density toward PZ Tel is 
not known, and there is insufficient information in 
the HETG spectrum to estimate it. For all the 
following analyses we have therefore assumed an 
hydrogen column density of lO^^cm"^, appropri- 
ate for the distance of PZ Tel (49.7 pc) and an av- 
erage interstellar medium density of 0.1 cm^'^. We 
have estimated that if the hydrogen column den- 
sity is as high as 10^°cm~^, the effect would be 
marginally observable only for the O and N lines 
at A > 18 A, with an absorption of about 5-10%. 

3.1. Emission Measure Distribution and 
Abundance Reconstruction 

We have derived the emission measure distribu- 
tion from individual line fluxes: one of the main 
advantages of this procedure, in comparison with 
a global fitting approach, is that we can use only 
the spectral information that we consider most re- 
liable. This choice involves the question of which 
lines to include in the em{T) reconstruction and 
which to exclude. In principle one would expect 
that lines which form over the same tempera- 
ture range should give consistent emission mea- 
sure values, and therefore different line selections 
should not provide significant differences in the 
final em{T). However, in some cases the above 
assumption will not be correct because of uncer- 
tainties in the atomic data, hidden line blending, 
and other line measurement difficulties. 

We have explored two diff'erent choices of spec- 
tral lines, which we describe below. For both 
of these selections we derived the emission mea- 
sure distribution using the method proposed by 
Kashyap & Drake (1998), which performs a search 
in the em(T) parameter space using a Monte Carlo 
Markov Chain method, with the aim of maximiz- 
ing the probability of obtaining equal model and 
observed line fiuxes. Moreover, we have also used 
two different methods for the abundance determi- 
nation so that the two approaches are effectively 
independent. One of the aims of this is to com- 
pare the results from the different approaches in 
order to find which results are the least affected 
by the methods and hence more robust. Hereafter 
we will refer to the first selection as Approach 1, 
and to the second as Approach 2. 
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3.1.1. Approach 1 

The first Approach is based on the use of aU 
lines strong enough to yield a useful flux measure- 
ment, and that are also reasonably isolated, den- 
sity independent and not affected by identification 
problems. With this choice the em{T) analysis 
is based on the most reliable line fluxes among 
those measured. Note that, for a solar abundance 
corona, this selection criterion will produce a set of 
measurements mainly composed of iron lines. For 
the MEG spectrum of PZ Tel the lines included in 
this selection are indicated in the last column of 
Table 1. 

With this choice we are able to derive, besides 
the em(T), the abundances, relative to that of 
iron, of all the chemical elements for which at least 
one line has been included in the selection. The 
abundance reconstruction is performed simultane- 
ously with the em{T) as follows: we first used only 
the selected iron lines to obtain an initial em{T); 
we then added lines of other elements, one element 
at a time, in order to obtain a new em{T) and the 
abundances of the involved elements relative to Fe. 
At each step we kept the previous set of ern(T) 
and abundances as starting point for the new re- 
construction. The sequence of elements has been 
estabilished in order to gradually extend the tem- 
perature range where the em(T) is constrained. 
Finally, the Fe abundance has been evaluated by 
comparing the observed continuum with the pre- 
dictions based on different metallicity values (see 
Argiroffi, Maggio, & Peres 2003, for more details). 

In principle, this selection criterion allows us 
to reach the finest temperature resolution in the 
em(r) because of the large number of lines in- 
volved. While many of the lines have emissiv- 
ity functions that peak at the same Tmax, these 
emissivity functions can have significantly differ- 
ent shapes and hence can provide more informa- 
tion on the structure of the em{T) than a single 
line at the same temperature. Of course, this as- 
sumption does depend on there not being signif- 
icant errors in the emissivity functions of the in- 
cluded lines that might instead lead to distortions 
of the derived em(T) relative to its true form. 

3.1.2. Approach 2 

The second Approach is based on the use of 
lines produced by He-like and H-like ions and 



has been proposed by Schmitt & Ness (2004) and 
J. J. Drake, in preparation. This choice is moti- 
vated by the simplicity of these ions and that they 
are likely to have more reliable emissivity func- 
tions than lines produced by ions with a larger 
number of electrons, and in particular the Fe 
lines. Note that by adopting this selection it be- 
comes necessary to separate the em(T) reconstruc- 
tion from the abundance evaluation because these 
quantities are more critically correlated than in 
the case of Approach 1. This is because there is no 
element, like iron in Approach 1, whose lines cover 
the whole temperature range. In order to separate 
the em{T) and abundance reconstruction, we have 
used the flux ratios of the He-like triplet resonance 
line to the H-like Lya line for each of the different 
elements for which these lines have been signifi- 
cantly detected. In addition to these line ratios, 
we have also included three prominent Fe lines for 
a better sampling of the temperature range. In 
the case of PZ Tel the lines selected are marked in 
the last column of Table 1 . 

In applying this method we also included mea- 
surements of the continuum flux at 2.4-3.4 A, 5.1- 
6.1 A, and 19-20 A. Fluxes were obtained by sum- 
ming the counts in each interval. We have checked 
that no strong line features were present in these 
intervals, in both the observed and predicted spec- 
tra. These continuum measurements serve two 
purposes. Firstly, they provide additional temper- 
ature information, based on the fall-off of the dom- 
inant bremsstrahlung continuum toward shorter 
wavelengths. Secondly, they provide a normaliza- 
tion for the em(T) since our approach based on 
line ratios described above provides only an esti- 
mate of the shape of the function. An absolute 
measurement of em{T) also enables us to express 
element abundances relative to hydrogen. Ele- 
ment abundances were determined when a satis- 
factory em,{T) had been derived by comparing pre- 
dicted and observed line fluxes for each element. 

This line ratio method was developed by 
J. J. Drake, in preparation, in order to investigate 
coronal element abundances in a large sample of 
stars. The H- and He-like resonance transitions of 
the abundant elements N, O, Ne, Mg and Si are the 
brightest lines in coronal spectra and so should be 
detected in all stars for which well-exposed spec- 
tra have been obtained. One advantage of this 
method is, then, that the same analysis using the 
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Fig. 3. — Emission measure distribution em{T) of 
PZ Tel derived with the two criteria introduced in 
§ 3.1. 

same hnes can be performed on many stars, hmit- 

ing any systematic bias caused by star-to-star line 
selection differences. One disadvantage, at least 
in principle, is that the line emissivities of He-like 
and H-like ions have in general a larger width in 
temperature than those of iron linos, for exam- 
ple. Any fine structure in the ern{T) function will 
therefore be more likely to be lost in this method 
than in an approach using a more extensive line 
list, such as our Approach 1. 

4. RESULTS 

4.1. Emission Measure Distribution and 
Abundances 

Considering the formation temperatures of the 
available spectral lines listed in Table 1, we have 
explored the em{T) in the logT range 6.0 — 7.5, 
with a bin size of AlogT = 0.1. 

In Fig. 3 we show the em{T) obtained with the 
two different methods and line lists, while in Ta- 
ble 2 and Fig. 4 we report the derived abundances. 
As explained in § 3.1, Approach 1 compares the 
observed and predicted continuum level in order to 
evaluate the iron abundance and hence the overall 
plasma metallicity. Thus we have not assigned an 
error bar to the Fe abundance in Fig. 4. This com- 
parison yields an iron abundance of ~ 0.3 times 
the currently accepted solar value (Grevesse & 
Sauval 1998). Note that the iron abundance so 
derived plays the role of a global factor both in 
the em(T) and in the abundance values derived 



Fig. 4. — Abundances in solar units (Grevesse & 
Sauval 1998) vs. elements sorted by increasing 
First Ionization Potential (FIP) for both the cri- 
teria. 

with Approach 1. 

The O, Ne, Mg, Al, Si, S, and Fe abundance es- 
timations in Approach 2 were performed using the 
lines that were used in the em{T) derivation; the 
N, Na, and Ni abundances were instead evaluated 
using their measured line counts listed in Table 1 
that were excluded from the em(T) analysis. 

The two emission measure distributions ob- 
tained with the two approaches (see § 3.1.1 and 
3.1.2) are compatible within statistical imccrtain- 
ties in almost all the bins of the explored tem- 
perature range. Both ern(T) distributions show 
a maximum at logT = 6.9, which is a common 
feature for active stars (e.g. Drake 1996; Sanz- 
Forcada, Brickhouse, & Dupree 2003a), and also a 
significant emission measinc for logT > 7. How- 
ever, some differences are present, mainly near 
logT ~ 6.5 and in the high temperature tails of 
the distributions. 

Several causes might be responsible for these 
differences. Errors in the atomic data may yield 
different solutions if different line sets are used. 
This is one of the primary motives of our Ap- 
proach 2, which should introduce minimal errors 
arising from atomic data. 

Another problem is due to the fact that the 
observed lino fluxes are related to the em{T) 
function through a Fredholm equation of first 
type, which has no unique solution (e.g. Craig 
& Brown 1976). Difi'erent em(T) and abundance 
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Table 2 

Element Abundances Ax = Nx /Nn . 



Approach 1 Approach 2 



Elem. 


Ax/AxQ 


Ax 




Ax/AxQ 


Ax 




Na 


1.2 ±0.4 


(2.6 ±0.9) X 10- 


-6 


1.5 ±0.4 


(3.3 ± 1.0) X 10- 


-6 


Al 


< 0.3 


< 1.0 X lO"'^ 




0.6 ±0.3 


(1.7 ±0.7) X 10- 


-6 


Ni 


< 0.5 


< 9.5 X 10"^ 




0.5 ±0.4 


(9 ±6) X 10-^ 




Mg 


0.28 ±0.03 


(1.1 ±0.1) X 10" 


-5 


0.34 ±0.03 


(1.3 ±0.1) X 10- 


-5 


Fe 


~ 0.3 


~ 9.7 X 10"*^ 




0.27 ±0.03 


(8.5 ±0.9) X 10- 


-6 


Si 


0.22 ±0.03 


(7.9 ±0.9) X 10- 


-6 


0.30 ±0.04 


(1.1 ±0.1) X 10- 


-5 


S 


0.27 ±0.10 


(6 ±2) X 10-6 




0.38 ±0.13 


(8 ±3) X 10-6 







0.47 ±0.06 


(3.1 ±0.4) X 10- 


-4 


0.41 ± 0.08 


(2.7 ±0.5) X 10- 


-4 


N 


< 1.1 


< 9.4 X 10"^ 




< 1.2 


< 9.7 X 10"^ 




Ne 


0.57 ±0.04 


(6.9 ±0.5) X 10- 


-5 


0.76 ± 0.06 


(9.1 ±0.8) X 10- 


-5 
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Fig. 5. — Comparison between observed line fliixes 
and those predicted by the ern(T) obtained with 
Approach 1. Wc indicate with the crosses the hne 
fluxes selected for Approach 1 and not for Ap- 
proach 2, while the filled circles represent the line 
fluxes included both in Approach 1 and 2. 

sets might then reproduce equally well the ob- 
served line fluxes within measurement uncertain- 
ties. This non-uniqueness might also be affected 
by the resolution and amplitude of the tempera- 
ture grid assumed a priori for the em{T) recon- 
struction. 

The two different line choices and methods pro- 
vide consistent abundance values in most cases, 
indicating that the derived abundances do not de- 
pend on small differences in the emission measure 



Fig. 6. — Comparison between observed line fluxes 
and those predicted by the em(T) obtained with 
Approach 2. The crosses and flUed circles are de- 
fined as in Fig 5. 

distributions. We find a significant discrepancy 
between the two models only for Si and Ne abun- 
dances. These small abundance differences are 
probably due to the slightly different values of the 
em{T) obtained with the two approaches in the 
temperature range where Si and Ne lines are pro- 
duced. 

In Fig. 5 and 6 we show the comparison be- 
tween observed and predicted line fluxes for both 
the solutions. From these plots we deduce that, 
as expected, the em{T) derived with Approach 1 
shows a better overall description of the observed 
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Table 3 

Electron Densities Calculated from the He-like Triplets. 



Ion 




//' 


log AV (cm--'') 


Ne IX 


6.6 


3.0 ±1.4 


< 11.8 


Mg XI 


6.8 


1.7 ±0.8 


12.6 ±0.6 


Si XIII 


7.0 


3.5 ±2.1 


< 13.5 



line fluxes of PZ Tel than in the other case, but the 
predicted line flux values of the resonance line of 
the He-like ions and Lya line of the H-like ions 
are more discrepant than those computed with 
Approach 2 solution (filled circles in Fig. 5 and 
6). Moreover, we note that the Approach 2 so- 
lution describes spectral features formed at lower 
temperatures somewhat better, as represented by 
a smaller spread between observed and predicted 
fluxes for logT < 7.0. The Approach 1 solution, 
on the other hand, shows a smaller spread between 
observed and predicted fluxes for logT > 7.0. 

Figure 4 shows that the PZ Tel coronal abun- 
dances, relative to solar photospheric values 
(Grevesse & Sauval 1998), do not show a strong 
FIP dependence, but rather show a pattern vs. 
FIP similar to that of other active stars (e.g. 
Brinkman et al. 2001; Drake et al. 2001; Audard 
et al. 2003; Sanz-Forcada et al. 2003b): elements 
at low and high FIP appear enhanced with respect 
to the medium FIP elements. It is worth noting 
that we cannot compute abundance ratios relative 
to photospheric values because those of PZ Tel 
are not known; therefore the real abundance vs. 
FIP pattern may differ significantly from the one 
plotted in Fig. 4 (e.g. Sanz-Forcada, Favata, & 
Micela 2004). 

Finally we note that, using the flux measured 
for the Na XI Lya line at 10.02A, we have per- 
formed — to our knowledge — the first coronal 
Na abundance estimate based on X-ray spectra. 
The Na abundance is useful for constraining the 
abundance vs. FIP pattern at low FIP. The iden- 
tification and subsequent flux estimation of these 
lines is nontrivial because this spectral region also 
contains a number of Fe XX and other lines. Ayres 
et al. (2001) and Phfllips et al. (2001) identified 
the Na XI Lya line in the HETGS spectra of 
HR 1099 and Capella, respectively. Sanz-Forcada 



et al. (2003b, 2004), in the HETGS of AB Dor 
and A And respectively, measured the flux of a 
spectral feature at 10.02 A, but they did not iden- 
tify it because of a lack of Na lines in the As- 
trophysical Plasma Emission Database (APED). 
Moreover, in the HETGS spectrum of AR Lac, 
analyzed by Huenemoerder et al. (2003) using the 
APED database, this emission line at 10.02 A is 
clearly visible and does not appear to be easily 
attributable to other species, such as Fe XX. In 
conclusion, we suggest that the spectral feature at 
10.02 A observed in several coronal X-ray spectra 
is to be identified with the Na XI Lya line, listed 
in the CHIANTI database. 

4.2. Electron Density 

Estimates of the electron density can be ob- 
tained using He-like triplets (Gabriel & Jordan 
1969). These spectral features consist of three 
lines (resonance r, intercombination i and forbid- 
den /) produced by transitions from n = 2 to 
n = 1 levels in He- like ions. For the physical con- 
ditions of coronal plasmas, the ratio f/i is pre- 
dominantly sensitive to density, while the ratio 
{i + f)/r is predominantly sensitive to tempera- 
ture. Note that the density and the temperature 
inferred from these lines will be weighted averages 
of for region in which the relevant triplet is formed. 

We have identified and measured in the MEG 
spectrum of PZ Tel the He-like triplets of Ne IX, 
Mg XI and Si XIII ions, while for the O VII, Al XII 
and S XV ions only the resonance line is measur- 
able (see Tabic 1). In Table 3 we report the f/i 
ratios, together with the derived electron densities 
Nc and la uncertainties. Note that, because of 
the unresolved blends of Ne IX i and / lines with 
Fe XIX lines, we have evaluated the Ne IX f/i ratio 
by subtracting the predicted fluxes of the Fe XIX 
lines from the measured values. We have obtained 
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density values from the Mg XI, while from Ne IX 
and Si XIII we have derived only an upper limit. 
We also point out here that our measured //i ra- 
tios have uncertainties of more than 40% and are 
all consistent with their respective low density lim- 
its at the 2<T level. 

5. DISCUSSION AND CONCLUSIONS 

It is interesting to compare the coronal proper- 
ties of PZ Tel found here with those of TW Hya 
and AB Dor, which are slightly younger and older, 
respectively, than PZ Tel. TW Hya is a classi- 
cal T Tauri star with a circumstellar disk and an 
age of about lOMyr (Webb et al. 1999). AB Dor 
is a Kl single star that has almost arrived on 
the main sequence. Collier Cameron & Foing 
(1997) have estimated the age of AB Dor to be 
in the range from 20 to 30 Myr, while Favata 
et al. (1998) derived a slightly older age, about 
> 35 Myr. In Fig. 7 we show the positions on 
the HR diagram of these three stars and have also 
superimposed the evolutionary tracks from Ven- 
tura et al. (1998) with Z = 0.01. The values of 
My and B — V have been derived from the Hip- 
parcos catalog (Ferryman et al. 1997). Note how- 
ever that TW Hya cannot be placed accurately on 
the HR diagram because of the large variability 
observed in its color indices (Rucinski & Kraut- 
ter 1983; Mckkadcn 1998). Rucinski & Krautter 
(1983) found that the spectrum of TW Hya is com- 
patible with a K7V star, and we have therefore as- 
sumed a B — V color index of 1.30 (Johnson 1966). 
The positions of these stars on the HR diagram 
(Fig. 7) confirm that TW Hya and AB Dor are in 
evolutionary stages just preceding and following 
the stage of PZ Tel. 

The coronal emission of TW Hya and AB Dor 
has been studied via Chandra HETGS observa- 
tions by Kastner et al. (2002), Sanz-Forcada ct al. 
(2003b) and D. Garcfa-Alvarez, in preparation. 
Comparison with our present results shows that 
PZ Tel is very similar to AB Dor in its emis- 
sion measure distribution, element abundances 
and electron densities. TW Hya, however, appears 
to be very different in many respects. Its emission 
measure distribution peaks at logT = 6.5, and 
there is no evidence of plasma at logT > 7.0 in 
the observed spectrum, while PZ Tc^l and AB Dor 
have their maxima at logT ~ 6.9, and both ex- 
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Fig. 7.— HR diagram for PZ Tel, TW Hya and 
AB Dor. Solid lines represent the evolutionary 
tracks of Ventura et al. (1998) evaluated for Z = 
0.01. Dotted lines represent the isochrones at 10^, 
lO'^ and 10^ yr. 

hibit considerable emission measure at higher tem- 
peratures. There are also abundance differences, 
especially regarding the Ne/Fe ratio. In TW Hya, 
Ne/Fe is larger than the solar value by a factor 
~ 10, whereas for AB Dor and PZ Tel the en- 
hancements are less extreme, with factors of ~ 5 
and ~ 2, respectively^. 

Finally, the electron densities of the emitting 
plasma of TW Hya, derived from the O VII and 
Ne IX He-like triplets, turn out to be much higher 
than those of PZ Tel and AB Dor. The differences 
in densities between PZ Tel, AB Dor and TW Hya 
can be appreciated in Fig. 8, where we show their 
observed MEG spectra in the region of the Ne IX 
lines. We show the Ne IX triplet region because 
it is the only triplet for which the observed spec- 
tra have sufficient S/N for all three stars (in fact, 
the O VII intercombination and forbidden lines in 
the spectrum of PZ Tel are not measurable). It 
is known that the Ne IX triplet is complicated by 
blends with highly ionized Fe lines that are clearly 
visible in the PZ Tel spectrum, e.g., at 13.51 A (see 
also Table 1 for more detail). The Fe XIX line at 
^ 13.51 A in the TW Hya spectrum is instead not 
observable because of the peculiar Ne/Fe abun- 
dance ratio and the relatively low temperature of 

^For the Ne/Fe abundance ratio of PZ Tel, we have used re- 
sults obtained with Aijproach 1. Approach 2 gives a Ne/Fe 
abundance ratio ~ 3, which also tends to confirm the ob- 
served trend. 
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Fig. 8. — Comparison of the Ne IX spectral region 
of PZ Tel, AB Dor and TW Hya. For each star 
we have plotted the first order spectra of Chandra 
MEG. 

the peak in em{T) for this star. However, it is still 
quite clear that the relative intensities of the i and 
/ lines for PZ Tel and AB Dor are very different 
from those of TW Hya. In fact the f /i ratio is 
2.9 ± 0.3 for AB Dor, 3.0 ± 1.4 for PZ Tel and 
0.44 ±0.12 for TW Hya, implying electron density 
values of logA^e ~ H, < 12, and ~ 13, respec- 
tively. 

As noted in § 1, Kastncr et al. (2002) attributed 
the dominance of plasma at a temperature of 
3 X 10^ K, with high electron densities of order 
log A'o '--^ 13 in TW Hya, to accretion rather than a 
magnetically heated corona. In the accretion sce- 
nario, the plasma is heated in a shock at the stellar 
surface. The strong difference wc sec here between 
the X-ray emission of the accreting TW Hya and 
non-accreting PZ Tel supports these conclusions. 
However, there are two outstanding questions that 
the accretion scenario does not answer: (i) why is 
the Ne/Fe abundance in the shocked plasma sim- 
ilar to that seen in some coronally active RS CVn 
systems (see e.g. Drake et al. 2001; Huenemo- 



erder, Canizarcs, & Schulz 2001) and so different 
from the presumed value in the accreting circum- 
stellar material? And (ii) why is there very little 
trace of underlying coronal activity similar to that 
in PZ Tel? 

If the Ne/Fe abundance seen in TW Hya is 
indeed representative of the ac;creting material, 
then significant compositional fractionation must 
be occurring in the circumstellar disk or in the 
inner disk-magnetosphcrc region from which the 
accreting plasma falls onto the star. The lack or 
weakness of "normal" coronal activity on TW Hya 
might logically be associated with accretion ac- 
tivity or the presence of substantial disk material 
close in to the star itself. In this respect, however, 
it seems that TW Hya might bo urmsual among ac- 
creting CTTSs. X-ray studies at low resolution of 
large samples of PMS stars (Tsujimoto et al. 2002; 
Nakajima et al. 2003) have indicated that CTTSs 
coronal spectra tend to be dominated by plasma 
at logT ~ 7.5 — much hotter than in the case 
of TW Hya. Note however that CTTSs are often 
highly absorbed and therefore the high difference 
in temperature could be the result of a selection 
bias because of their distances larger than that of 
TW Hya. 

Our study of PZ Tel shows it to be much more 
similar in ern(T), abundances and electron den- 
sities to AB Dor than to TW Hya. Because 
PZ Tel is younger than AB Dor, our results in- 
dicate that coronal properties do not change ap- 
preciably during the evolution of PMS stars from 
the post T Tauri phase to the main sequence. On 
the other hand, if TW Hya is assumed to be rep- 
resentative of CTTSs, it follows that as soon as 
disk dissipation and significant accretion processes 
have ended, the coronal structure quickly evolves 
to that of any other star of similar spectral type 
and rotation rate. 
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